We suggest that high-order g modes can be used as a probe 
of the internal magnetic field of SPB (slowly pulsating B) stars. 
The idea is based on earlier work by the authors hich analyti- 
cally investigated the effect of a vertical magnetic field on p and 
g modes in a plane-parallel isothermal stratified atmosphere. It 
was found that even a weak field can significantly shift the g- 
mode frequencies — the eff'ect increases with mode order In 
the present study we adopt the classical perturbative approach 
to estimate the internal field of a 4 solar mass SPB star by look- 
ing at its effect on a low-degree (I - I) and high-order (n = 20) 
g mode with a period of about 1.5 d. We find that a polar field 
strength of about 110 kG on the edge of the convective core is 
required to produce a frequency shift of 1%. Frequency split- 
Jjf^ • tings of that order have been observed in several SPB variables, 
in some cases clearly too small to be ascribed to rotation. We 
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suggest that they may be due to a poloidal field with a strength 



> 

o 



> 

(N 



in 
o 

o 



of order 100 kG, buried in the deep interior of the star 
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Abstract. 

1. Introduction 

It is well known that a magnetic field produces a fre- 
quency splitting of stellar oscillations (Ledoux & Simon 1957). 
This effect has been extensively studied, using a perturba- 
tive approach (e.g. Goossens 1972; 1976a,b; 1977; Goossens, 
Smeyers & Dennis 1976; and more recently by Gough & 
Taylor 1984; Dziembowski & Goode 1984, 1985; Gough & 
Thompson 1990; Shibahashi & Takata 1993). This method 
must be handled with care near the surface, where the magnetic 
pressure dominates over the gas pressure, and where the acous- 
tic modes therefore strongly couple with the Alfven modes 
(e.g. Biront et al. 1982, Roberts & Soward 1983, Campbell & 
Papaloizou 1986, Dziembowski & Goode 1996); such regions 
are better treated using a non perturbative treatment (Bigot et 
al. 2000). 

In an earlier paper Hasan and Christensen-Dalsgaard 
(1992) analytically detemuned the frequency shift of p and g 
modes in an isothermal plasma due to a homogeneous verti- 
cal magnetic field. Using the full MHD equations, they found 
that even a weak field (more precisely when » 1, where /3 
is the ratio of gas to magnetic pressure) can produce a signifi- 
cant shift of g-mode frequencies, while the effect on the p-mode 
spectrum is comparatively small. In principle this means that g- 
mode frequencies offer a diagnostic to probe the internal field 
of stars in which g modes have been observed on the stellar 
surface. 

Extensive observation campaigns have uncovered the ex- 
istence of a class of variable stars known as slowly pulsating 
B (SPB) stars which are multiperiodic typically over a time 
scale of days (Waelkens 1991; De Cat et al. 2005 and refer- 
ences therein). These pulsations have been identified with low 
degree / (typically / = 1 and 2) g modes of high order, that are 
excited by the k mechanism in the metal opacity bump at a tem- 
perature of about 2 X 10^ K (Dziembowski et al. 1993). These 
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modes often occur in multiplets with closely spaced periods 
(with a typical separation of 1%). In some cases this separation 
can clearly not be due to rotational splitting, which would yield 
much larger spacings, as was pointed out by De Cat and Aerts 
(2002). In this letter we propose that such frequency splittings 
are due to the presence of a magnetic field. If this hypothesis 
is correct, then the splitting of frequencies can be used to esti- 
mate the field strength in the interior of SPB stars. 

2. Magnetic frequency splitting of g modes 

As was established in the early papers quoted above, the fre- 
quency shifts 6a) due to a magnetic field are given by: 

6oj 1 r-[(VxB')xBH-(VxB)xB']-^W 
— = ■ (1) 

here w is the angular frequency, ^ is the Lagrangian displace- 
ment with radial and horizontal components and ^h, respec- 
tively, and B' = Va(^AB) is the perturbation in the equilibrium 
magnetic field B. We shall deal here with fields whose energy is 
sufliciently small compared to the gravitational energy, so that 
we may neglect the structural changes caused by the Lorentz 
force. 

Let us consider a poloidal axisymmetric field of the form: 



B=Bn 



2b{x) cos e, --dA^b) sin 6, 



(2) 



where x = r/R is the normahzed radial coordinate, R is the 
radius of the star and bix) = 0(1), so that Bq characterizes the 
field strength. The displacement ^ can be expressed as: 



ITfl 

sm0 



exp iojt , (3) 



where Y" is the spherical harmonic of degree I and azimuthal 
order m (henceforth, denoted as F). For a poloidal field given 
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by equation (2), the perturbed field B' is: 



[/(/+ l)cos0y H-sine^eF] 



1 .,...2n..^r. 



5fl(sin^ 9 Y) ^dAx^b) 
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cosOdgY -d,(2x^hb) - sin0Y-dJ^rdAx^b)] 

X X ^ ' 



.dAx'b) 



, ^0 



cos 6 1 



sin0 X 



smU X- 

Y-dA2x^i,b)-deY%dA^b) 
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Table 1. Frequency and periods for g modes of different radial 
orders and degree Z = 1 , for a 4 star with an age of 94 Myr. 
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0.788 
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1.066 
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29.31 


0.395 


10 


25.15 


0.460 



For g modes of high radial order n, it is straightforward to show 
that in equation (1), the dominant term in the integrand of the 
numerator is the first one: 



- [(V X B') X B] ■ f = 



(4) 



|B'|- 



2 d 

X ax 



cos 9 



80 



+ m' 



cos 6 ^ 
Y 



sin0 



whereas in the denominator it is |^hP, since » for n » 1. 
In integrating by parts, we have neglected the surface terms; 
this region anyhow requires a special, non-adiabatic and non- 
perturbative treatment (see for instance Bigot et al. 2000). From 
equation (1) it therefore follows that 
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where pc is the central mass density, 
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Obviously these constants do not depend on the sign of m: the 
magnetic field reduces the (2/ + 1) degeneracy of the eigen- 
modes to only /+ 1, as already pointed out by Ledoux & Simon 
(1957). For / = 1 and 2, we find: 



Cu. - 3 



C 



1,1 



C2 2 — C2 -2 — TT 

21 



It is convenient to express equation (|5j as: 



— - ScBq, where Sc - 



Clm I 



(8) 
(9) 

(10) 



which we henceforth refer to as the splitting coefficient. This 
coefficient increases rapidly with period, since I increases also. 



Note that a toroidal field would produce a much lesser split- 
ting, since the leading term in equation would then be |^hP 
instead of |5t^hP- 

Finally, let us recall that the rotational splitting for g modes 
of high order is given by 



1 



1 



/(/ + 1) 



Q. 



(11) 



where we have again assumed that \^r\ ^ l^hl; we see that the 
frequency spacing is of the order of the average angular veloc- 
ity Q. 

3. Model 

Following Dziembowski et al. (1993), we consider a model 
SPB star with the following parameters: M - 4Mq, 
logiL/Lo) = 2.51, log Teif = 4.142, = 0.37. An equilibrium 
model for such a star was calculated using the Aarhus Stellar 
Evolution Code (ASTEC) (e.g. Christensen-Dalsgaard 1982, 
1993); this used the Eggleton et al. (1973) equation of state and 
OPAL opacities (Iglesias & Rogers 1996) and ignored diffusion 
and settling. The evolved star had an age of 94 Myr We calcu- 
lated the g modes of the above star using the Aarhus adiabatic 
oscillation package (e.g. Christensen-Dalsgaard & Berthomieu 
1991). Table 1 lists the cyclic frequencies v and periods P for 
high-order g modes corresponding to Z = 1 and 2. 

We first evaluate the frequency shift due to a magnetic field 
for the mode (i.e. a g mode of radial order 20 and degree 
I - 1), which has a period of 1.5 d (a typical period for a SPB 
star). The horizontal eigenfunction (^h) of this mode is shown 
in Figure 1 : the radial component is normalized to unity on 
the surface r = R of the star 

We calculated the numerator of equation (|6j with different 
functional forms b{x) for the magnetic field. When the field 
is constant throughout the star, the main contribution to that 
integral comes from above x w 0.8, where has its largest 
amplitude. But when the field is buried below that depth, or 
when it tapers off as oc x"* with q > 1, its main contribution 
originates from a small region just above the convective core 
at Xc - 0.106, where there is a steep gradient in the molecular 
weight in this evolved star. When we choose bix) - (x/xcy. 
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Fig. 1. Horizontal component of the eigenfunction (scaled 
such that ^j. = 1 at the surface) as a function of normalized ra- 
dius {rlR) for a g^^ mode with a period of 1 .5 d for a 4 M© SPB 
star of 94 Myr The inset shows an enlargement of the region 
close to the edge of the convective core, where the eigenfunc- 
tion displays strong oscillations due to the presence of a steep 
composition gradient. 

Table 2. Splitting constant (5 c) and value of the polar field 
(Bpoi), at the edge of the convective core for / = 1 g modes 
that would produce a frequency splitting of 1 %, for a 4 Mq star 
with an age of 94 Myr. 



Mode 


P (day) 


(Gauss-2) 


Spoi(kG) 




1.497 


2.278 X lO^*^ 


111 


sis 


1.159 


3.765 X 10-' 


274 




0.789 


2.321 X 10-** 


1102 



the results depend little on q; those presented hereafter were 
obtained with q -!>. 

For this mode a? - 2.36 x 10"' s"^. Using Ci,i - Ci,o = 
1/5, one finds 5c = 2.278 x 10"^. From equation (O, one 
deduces that in order to produce a 1 % frequency shift in 6(jj / w, 
the polar field just above the convective core has to be Bpoi - 
llOkG. 

For comparison, we consider now the effect of the field on 
a g mode of order 10 and the same degree I - \. The hori- 
zontal component of the eigenfunction for this mode is shown 
in Figure 2. In this case, Bpoi - 1100 kG, which is an order 
or magnitude larger than for the n = 20 mode with the same 
degree. 

Tables 2 and 3 give the splitting coefficient 5 c for / = 1 and 
I - 2 g modes of various orders, and the corresponding polar 
field strength Bpoi required to produce a 1% frequency shift. 
In Table 3 we separate the contributions due to C20 - C22 and 
C21 - C22 terms which produce different frequency shifts. Note 
how rapidly Sc increases with radial order n. 



Table 3. Splitting coefficient (Sc) and value of the polar field 
(Bpoi), at the edge of the convective core for I - 2 g modes 
corresponding to different transitions that would produce a fre- 
quency splitting of 1%, for a 4 Mq star with an age of 94 Myr. 



C2.0 - C2.121 £"2.11 1 - C2.121 

Mode P(day) 5c (G'^) Bpoi(kG) 5^ (G'^) Bpoi(kG) 

1.497 6.435 xlO^' 210 4.827 x 10^' 242 

g5o 1.281 3.422x 10-' 288 2.566x 10-' 332 

A 1.066 7.897x 10-' 1890 5.922x 10-' 2180 
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Fig. 2. Same as in Fig. 1, but for the ^[q mode. The inset shows 
an enlargement of the region close to the convective core. 

4. Discussion 

Our calculations suggest that high-order g modes can provide a 
sensitive diagnostic of the internal magnetic field in SPB stars, 
at least when these are sufficiently evolved. The main contribu- 
tion to the splitting coefficient then comes from a small region 
just above the convective core (close to x = 0.106 in a 4 
star of 94 Myr), where a steep helium gradient has built up 
due to the shrinking of the convective core. This causes a sharp 
peak in the Brunt-Vaisala frequency A^, which is reflected in a 
series of closely spaced nodes in the horizontal eigenfunction. 
In order to examine the impact of such a region, we repeated 
the calculation for a ZAMS star with the same chemical com- 
position and physical parameters as our evolved SPB star, and 
found indeed that the frequency splitting was then much less 
sensitive to the strength of the magnetic field. 

We concentrated on g modes of low degree and high order, 
which are typically excited in SPB stars. We found that a polar 
field of 1 10 kG in the vicinity of the convective core causes a 
splitting of 1% for a gj^ mode. For a simple dipolar configu- 
ration of the magnetic field, this would translate into a 120 G 
polar field on the surface of the star; we quote this figure only 
for illustration purpose, since there is no way to deduce the sur- 
face field from the deep field. We have checked that the rest of 
the star contributes httle to the splitting, provided that the field 
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is buried below the depth of about x = 0.80, or that it tapers 
off at a faster rate than b ocl/x.ln such a situation we find that 
high-order g-modes can be be used to probe the deep interior 
field. This result is not sensitive to the precise field configura- 
tion. 

With moderate- to low-order modes, the diagnostic is much 
less sensitive. For a I - 2, n = lOg mode with a period of 
1 . 1 d, a polar field of about 2 MG would is required to produce 
a frequency splitting of 1%. 

So far only one SPB star has been detected with a magnetic 
field. Recently, Neiner et al. (2003) have reported the discovery 
of a field on the SPB star f Cas, in which a non-radial pulsation 
with a period P - 1.56 d was detected. A field strength for 
the time-averaged line of sight polar component of 330^lf G 
was inferred. If this field is the visible part of a deeply rooted 
magnetic field, it could also leave a signature in the splitting 
of high-order g modes. We should emphasize that according to 
our results, a field of order 100 kG at the edge of the convective 
core would be required to produce a 1 % splitting in typical g 
modes. 

We considered here only a purely poloidal configuration, 
similar to many earlier papers quoted in Sect. 1. Such a config- 
uration is known to be unstable. As was shown by Tayler and 
collaborators (Tayler 1973, Pitts & Tayler 1985), and as was il- 
lustrated recently by the numerical simulations of Braithwaite 
and Spruit (2004), the configurations which are likely to re- 
sist non-axisymmetric MHD instabilities are combinations of 
large-scale toroidal and poloidal fields of about equal strength. 
Taking this into account is unUkely to alter our conclusions 
concerning the detectability of the deep magnetic field, because 
the frequency splitting would be much more sensitive to the 
poloidal than to the toroidal component in such a combined 
field. 
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